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Abstract 

The BFKL and the unified angular-ordered equations are solved to determine the gluon 
distribution at small x. The impact of kinematic constraints is investigated. Predictions 
are made for observables sensitive to the gluon at small x. In particular comparison is 
made with measurements at the HERA electron-proton collider of the proton structure 
function ^(x, Q 2 ) as a function of InQ 2 , the charm component, F£(x, Q 2 ) and diffractive 
J/ip photoproduction. 
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1. Introduction 

Understanding the small x behaviour of the gluon density is one of the most challenging 
problems of perturbative QCD. It has become topical with the commissioning of the electron- 
proton collider HERA at DESY. Indeed experiments at HERA are probing the gluon density in 
the region x < 10~ 3 by observation of the behaviour of the proton structure functionP] F2 as a 
function of In Q 2 and by the measurements of diffractive J/ ip photoproduction, dijet production 
etc. 

First we recall that for moderate values of x, say x > 0.05, observable quantities such as 
F 2 (x, Q 2 ) are determined in perturbative QCD by the mass factorization theorem in which the 
collinear logarithmic singularities, that arise from gluon emissions in the partonic subprocesses, 
are absorbed into universal parton densities. The absorption of the collinear singularities make 
the densities "run" with a Q 2 dependence determined by the Altarelli-Parisi (GLAP) evolution 
equations, although the absolute values of the densities are not calculable in perturbative QCD 
but have to be input at some scale, say Qq. In fact Altarelli-Parisi evolution resums the lead- 
ing as \n(Q 2 /Ql) contributions where, in a physical gauge, the a s , ha n (Q 2 /Ql) contribution is 
associated with a space-like chain of n gluon emissions in which the successive gluon transverse 
momenta are strongly ordered along the chain, that is q\ x -C . . . -C q^ n <C Q 2 . The coef- 
ficient (splitting) functions and anomalous dimensions of Altarelli-Parisi evolution have been 
calculated to next-to-leading order. This corresponds to the situation in which a pair of glu- 
ons are emitted without strong ordering (and iterations of this configuration). Then the 
contributions contain a power of as unaccompanied by \ia(Q 2 /Ql). 

At sufficiently high electron-proton cm. energy, y/s, we encounter a second large variable, 
1/x ~ s/Q 2 , and we must resum the leading ashi(l/x) contributions. In this regime the 
dominant parton is the gluon and the key ingredients to calculate hard scattering observables 
are the /^-factorization theorem and the BFKL equation for the gluon distribution F(x, k^) 
unintegrated over its transverse momentum fc^ @, Hi- The BFKL equation, which sums the 
leading as ln(l/a;) contributions, may be written in the form || |7j 



F(x, k 2 T ) 



1 dz f d 2 qx 



F^(x,k 2 T )+a s [ - f 
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(1) 



with a~s = 3a,s / tt, and where F^ is the non-perturbative driving term of the equation. To gain 
insight into the equation we refer to Fig. 1 which shows the basic unit that, on iterating (|IJ), 
leads to a chain of sequential gluon emissions. These real emissions come from the first term 
under the integral, where kx + qr = k' T of Fig. 1. The second term corresponds to the virtual 
contributions which arise from gluon loops along the chain. The real and virtual contributions 

3 In deep inelastic scattering we do not probe the gluon directly, but instead via the process 73 — > qq. The 
longitudinal fraction x of the proton's momentum that is carried by the gluon is therefore sampled over an 
interval bounded below by the Bjorken x variable Xb = Q 2 /2p.q, where as usual p and q are the four momenta 



of the incoming proton and virtual photon respectively and Q 2 



-q 2 . We omit the subscript b since the 



distinction between x (of the gluon distribution) and Xb in F 2 (xb,Q 2 ) is evident. 
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cancel as qr — > to ensure that there is no singularity in (|l]) . The strong ordering in is now 
no longer applicable and instead we have a "random walk" or diffusion in In q\ as we proceed 
along the chain to gluons with smaller and smaller values of x. The enlarged phase space 
leads to a x~ x type growth of F as x —>■ 0. In fact, since the kernel of the BFKL equation is 
scale invariant, the solution can be found by Mellin transform techniques and given analytically 
for x — > (for fixed coupling as). In particular we have the famous result, obtained in the 
original BFKL papers ||, that 

A = a 5 41n2 (2) 

corresponding to the maximum of a continuous spectrum of eigenvalues of the kernel of the 
(Mellin-transformed) BFKL equation. 

The BFKL equation sums the leading order asln(l/a;) contributions. The next-to-leading 
terms are not yet known. However, there is one kinematical constraint which should be imple- 
mented, namely f?|, || 

k 2 -. > zq\. (3) 

The constraint^ arises since, in the small z regime where the BFKL equation is valid, we require 
that the virtuality of the exchanged gluons arises mainly from the transverse, rather than the 
longitudinal, components of their momentum; that is 

\kf « k'l (4) 



There are also constraints from energy-momentum conservation ||. For example, for the 
gluon of transverse momentum k' T forming one of the links of the BFKL chain we require 
(q + k') 2 > 0, where q is the four momentum of the deep inelastic photon probe. That is 

k'l < ^ ~ W 2 (5) 

x 

at small Bjorken x. However (|3]) gives an implicit bound on which is usually much more 
restrictive than (|5|). This can be seen by considering a given value of k\. Then a high value of 
q\ implies an equally high value of and (|3]) becomes 

k' 2 < ^. (6) 

Since z > x, this bound is tighter than (|^) except for low values of Q 2 which are much less than 
k\. This is fortunate because it is possible to study the effects of imposing (|) on the solution 
of the BFKL equation. Indeed, the introduction of the upper limit of integration, q\ < k^/z, 
preserves the scale invariance of the BFKL equation and therefore allows the Mellin transform 
technique to be used to find the modified form of the solution as x — > 0, see Section 2. The 
modifications enter at next-to-leading order and above. They are found to reduce the exponent 

4 Throughout we call this a kinematic constraint although its origin is partly dynamical. We sketch the 
derivation of (||) in Section 3. 
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A of the x~ x behaviour from the BFKL value given in (0). Besides finding the analytic solution, 
we also solve the BFKL equation numerically and present results for the effective value of A 
at non-zero values of x so that we can study how the asymptotic (x — > 0) analytic result is 
approached. 

If as is allowed to run then the integrand in the BFKL equation is weighted more to the 
infrared end, and the effects of the kinematic constraint (|J) will become weaker. We investigate 
the effect in terms of the unified CCFM evolution equation |7], ID], TI\, which reduces to the 



BFKL equation in the small x limit. The CCFM equation is based on the coherent radiation 
of gluons, which leads to an angular ordering in their emissions [E|. The gluon distribution, 



which satisfies the CCFM equation, depends on an additional scale that is required to specify 
the maximum angle possible for gluon emission, which is given essentially by the scale of the 
probe. The gluon distribution therefore becomes a function of three variables F(x, fcf., Q 2 )- In 
particular, angular ordering introduces a constraint 



4<^(l-z) 2 (7) 



into the integrand of the equation for the gluon. At small z, this angular-ordering constraint 
q\ < Q 2 /z 2 appears much weaker than the kinematic constraint q\ < k\jz of @. Clearly 
for Q 2 > k\ we would expect no effect from angular ordering and the gluon distribution would 
become independent of Q 2 . In Section 3 we solve the CCFM equation numerically, with and 
without imposing the kinematic constraint (|3|), so as to study these effects. 

In Section 4 we use the results for the unintegrated gluon distribution, together with the 
fcT-factorization theorem, to study the effect of the kinematic constraint on observable small 
x processes. In particular we make predictions for the behaviour of F2(x, Q 2 ) as a function 
of InQ 2 at small x, for its charm component F|(x,Q 2 ), and for high energy diffractive J/ip 
photoproduction. Data are available for these three processes, all of which are sensitive to the 
gluon at small x ^, We present our conclusions in Section 5. 

2. Constraints on the gluon from the BFKL equation with fixed as 

In this section we investigate the effect of imposing the kinematic constraint q\ < k\jz on 
the solution of the BFKL equation, ([]]). In particular we wish to solve the equation for the 
unintegrated gluon distribution which incorporates the constraint, 



F(x,k 2 ) = F^\x,k 2 T 




x -Ak T + qT\ 2 ) - 0(4-4)^(f,4) 



(8) 

and to determine the value of the exponent A of the x~ x behaviour of the solution. The solution 
can most easily be obtained by rewriting the equation in terms of the gluon distribution in 
moment space 

F(w,k%.) = f dxx"- 1 F(x,k%). (9) 



o 



Then the BFKL equation becomes 



F{u, 4) 
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irq% 
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z 
2\ eV. . r_2 
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F(u, \k T + qr\ 



- 9 (A£ - qr) F(oj, k\ 



(10) 



where the first B-function imposes the kinematic constraint on the real emissions. As usual we 
solve the equation by taking the Mellin transform, and introduce the variable p conjugate to 



1 ' ^ " hi 



c+ioo 



d P (k 2 T yF(cu,p). 



F(w, k\ 

then the BFKL equation reduces to the algebraic relation 

F(u,p) = F®(u>,p) + ^K(u,p)F(co,p) 

00 

where the kernel 



(11) 



(12) 



K{u,p) 



d 2 g T 
■nq\ 



\k T + qr\ 



8(4-4) + (f)" 



6(g2 -k 2 T )\- 6(4 - q 2 ) 



(13) 

The expression {. . .} arises from the {. . .} in ([10|), up to a factor of 1/uj. The kernel is scale 
invariant, as becomes evident once we make the substitution q\ = k\u. Then (|i~3"D becomes 



1 dv r 

5 , — (l + 2v^cos0 + M )" {l + u"- p ) - 1 

Z7T JO U 1 
1 dll r 

— \F(-p,-p;l;u)(l + u»-o) - 

u 1 



(14) 



where F(—p, —p; 1; u) is the hypergeometric function, and where the p term arises from the 
integration over the interval 1 < u < oo when we substitute u — > 1/u. 

The solution of the BFKL equation for the double Mellin transform F(u,p) has the usual 
form 

f (( V,p) 



F(u,p) 



(15) 



1 - (as/u})K(u,pY 

except that now due to presence of the q\ < k\j z cut-off K depends on u, as well as p. Inverting 
the moments in the standard way we find that the unintegrated distribution is given by 

rc+ioc 



F(x, k 2 T ) 



1 

2tH 



dp{k 2 T ) p R{p)x- ul{p) (16) 

where we have performed the u integration, and picked up the "leading" pole of F(uj, p) at 
uj = (jj(p), with residue R(p). From (^) we see that the pole position uj(p) is a solution of the 
equation 

(17) 



uj(p) 



The exponent A of the x~ x behaviour, or in other words the pole position uJ(p) which determines 
the leading behaviour of ( Jl6| ) as x — > 0, is the maximum value of oo(p) as p varies along the 
contour c — ioo to c + ioo where c lies in the interval (— 1, uj). The value ul(p) is found from 



dK(u(p),p) 



dp 







which occurs at the saddle point p = p 



c. The resulting values for the intercept 
A = U(p) 



(18) 



(19) 



are shown as a function of as by the continuous curve in Fig.2. 

The curve is to be compared with the leading-order intercept, shown as a dashed line, which 
is obtained from the BFKL equation without the kinematic constraint imposed. In this case 
we have uo = in K(u, p) in (O) and so K is a function of p alone 



K(u = 0,p) = 2^(1) 



(20) 



where \l/ is the logarithmic derivative of the Euler T function : ^(z) = T'(z)/T(z). Along the 
contour of integration in (|1(J) the pole position u = asK reaches its maximum value at the 
"saddle" point p — c — — |, and so we recover the familiar leading-order BFKL intercept 



A = a s K {uj = 0, p 



a s 4 In 2. 



(21) 



It is easy to see that the imposition of the kinematic constraint in the BKFL equation introduces 
higher-order corrections to the value of the intercept. For instance we note that the pole in 



( |15| ) is a solution of 



to 



askiyj = 0, p) + a>su 



dk(cu,p) 



duo 



+ 



(22) 



ui=0 



If we were to keep only the next-to-leading term in our determination of the intercept, then we 
would have 

A = S<j41n2(l -4.15a<?), 

which is shown by the dotted curve denoted by (c) in Fig. 2. The importance of the higher-order 
corrections is evident, even for values of as = 3a,s/7r ~ 0.1. 

The above analytic procedure yields the "asymptotic" or leading x~ x behaviour of the gluon 
distribution F(x, k^) as x — > 0. We may solve the modified BFKL equation numerically to see 
how the asymptotic value of the intercept A is approached as x decreases, for different values of 
k\. To be precise we solve (||]) numerically using for the driving term F(°'(x, k^) a conventional 
distribution of the form 3(1 — a;) 5 iVexp(— k\jk^) with k = 1 GeV. In Fig. 3 we show the 
effective slope A e fr, obtained from the solution F(x, k^) via 



A, 



cff 



<91nF 
<91n(l/x) 



(23) 
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as a function of x for different values of for a fixed S5 = 0.2. We see that the approach 
to the asymptotic {x — > 0) value, A = 0.35, depends on the value of k\. Since A is the leading 
singularity we would expect an approach from below. However, for large k\ the x dependence 
of the BFKL diffusion pattern, 



exp 



-A 



\xi 2 {k 2 T /k T ) 
ln(l/x) 
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overrides this behaviour and causes A e g to approach the asymptotic value from above. Here k T 
is associated with the non-perturbative input distribution. 

2.1 Kinematic constraint in the folded form of the BFKL equation 

In Section 3 we study the effect of the kinematic constraint in the more realistic case when 
as is allowed to run. There we solve the evolution equation in "folded" form. That is the form 
in which all the virtual corrections and all the unresolved real gluon emissions are resummed. 
Unresolved emissions are those with q\ < fi 2 , where fi 2 specifies the resolution. Imposing the 
kinematic constraint q\ < k\jz on the folded BFKL equation is a little more involved than it 
was in the "unfolded" form (B) in which the real and virtual terms appear on an equal footing, 



that is to the same order in as- To see this we first rewrite the unfolded BFKL equation, (]10 ) 
in the form 

d 2 q T 
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(24) 



where we have simply divided the integration according to whether q\ is less or greater than 
and carried out the z integrations. The kinematic constraint is responsible for the factor 
( ) w in the second integral. Now we divide the first integral according to whether q\ < fi 2 or 
q\ > fi 2 so that we can resum all the virtual and unresolved real terms. Providing fi 2 is not 
too large, we can evaluate the q\ < [i 2 part explicitly: 

d 2 q T 



Os_ 

to 



7cq% 



F(uJ,\k T + q T \ 2 ) &(fi 2 
- -^F(u,k 2 ) f 

UJ Ju 



Cfr) - F{uj } k 2 T ) Q(l4-qi) 



- F(u, k. 

UJ 




(25) 



where 



uj = asHk^/fi 2 ). (26) 

The result (53) is the residual virtual correction to F(uj, k^) which remains after the cancellation 
of the unresolved real and virtual singularities. Using (^Sj) , we see that then fl21|) becomes 

d 2 q T 



F(u, k 2 T ) 




Q(k 2 T - 4) Q(q 2 T - n*) F(u, \k T + q T \ 2 ) 



uj + UJ 



Q{q^ - k 2 T ) Q{q^ - fi 2 ) F(u, \k T + q T 



(27) 



where = ojF°^ = wF + to). We invert the moments and transform back to x space, 
and find 

F(x, k%) 



F(°\x,k 2 T ) + 
_ f 1 dz f d 2 q T 
itq\ 



a s / — 

fx Z 



'k 



Q\-±-zj Q(q 2 T - //) A R (z, 14, q 2 T , fi) F [—,\kr + q T 
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(28) 



where A R , often called the non-Sudakov form factor, is given by 
K R (z,kl,ql^ 2 ) 



cxp 



-as 



z> U \q 2 T 




?4 Q(k 2 - q' 2 ) 0( g > 2 - S 
q T 



(29) 



if q\ < h\ 



(zqfr/kr) u if <?t > k Ti 
where UJ is given by Q2TSD . The form factor screens the 1/z singularity in 

We see that the inclusion of the kinematic constraint q\ < k\jz in the real gluon emission 
part of the unfolded BFKL equation, (§), has the additional effect of modifying the non-Sudakov 
form factor in the folded equation (|28|) . The modification occurs for q\ > k\ and makes the 
form factor larger, but always, of course, satisfying A R < 1. 



3. Constraints on the gluon from the CCFM equation 



The BFKL equation, which resums the asln(l/x) contributions, is applicable at small x 
and moderate Q 2 , whereas at larger x and large (Q 2 /Ql) Altarelli-Parisi (GLAP) evolution, 
which resums the asln{Q 2 /Qq) terms, is appropriate. A theoretical framework which gives 
a unified treatment throughout the x, Q 2 region has been formulated by Ciafaloni, Catani, 
Fiorani and Marchesini JTIJ . The CCFM approach is based on the coherent radiation of gluons 
which implies angular ordering of the gluon emissions along the chain. The CCFM equation 
embodies both the BFKL and GLAP equations in the appropriate kinematic regimes. In the 
small x region the CCFM equation may be approximated by |l 1 1 



F(x, k 2 , Q 2 



F^(x,k 2 T ,Q 2 ) + 

" l dz r d 2 q 
z J nq 2 



0(Q - zq) A R {z, k 2 T , q*) F - \k T + q T \ 2 , q 



(30) 



where it is convenient to impose the angular ordering in terms of rescaled transverse momenta 

q T 



q 



1 - z' 



(31) 
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The angular-ordering constraint then becomes q < Q/z. Here the non-Sudakov form factor A# 
is given by 



A R (z,k 2 T ,q 2 ) = exp [ ^ / ^ 9(4 - g' 2 ) 9(g'-z'g)j . (32) 



Finally, note that we have allowed a$ to run in (|30|) . The solutions fcf,, Q 2 ) of this equation 
were studied in ref. [0. As in ref. fill] we restrict the transverse momenta of the gluons along 



the chain to be above 1 GeV 2 . As we shall see it turns out that this physically reasonable 
choice of infrared cut-off gives a satisfactory normalization of all the observables sensitive to 
the small x behaviour of the gluon distribution, F, with the exception of F£ at low Q 2 . 

The novel feature of the CCFM, as compared to the BFKL equation, is that the solution 
is dependent on Q 2 . The origin of the dependence comes entirely from angular ordering, 
q < Q/z. In the BFKL leading ln(l/x) limit F becomes independent of Q 2 and indeed the 
CCFM solutions exhibit this behaviour for large Q 2 . However, for small Q 2 non-leading ln(l/x) 
effects become important via the angular ordering constraint and F decreases with decreasing 

g 2 0. 



Here we study the impact of the kinematic constraint q\ < k\jz. To be precise the con- 
straint is actually 

q 2 T < (l-z)k 2 T /z. (33) 

We sketch the derivation. The key observation is that the virtuality k 2 of a gluon along the 
chain should arise mainly from the transverse, rather than the longitudinal, components of the 
momentum for the small x approximation to be valid. Now in terms of the light-cone variables 
k = k ± k 3 

k 2 = k + k- - k 2 T (34) 

so we require 

k 2 T > \k + k~\. (35) 

From Fig. 1 we see that 

k~ = k'~ - q~ ~ -q~ = ~ql/q + , (36) 
where the last equality follows from the on-shell condition for the emitted gluonQ Thus 

k+k " ^ -£ = -*£f = ~~r~z * (37) 

The kinematic constraint ( |33] ) then follows directly from fl3"7p and (^). 



In ( |3q ) and (37) q are the light-cone components of the 4- momentum of the emitted gluon, whereas 
elsewhere in this section q denotes the rescaled transverse momentum defined by (|3l|). 



8 



Here we wish to compare the gluon distribution F(x, k^,Q 2 ) obtained by solving (|30D , 
with the solution obtained if the kinematic constraint ( |33| ) is imposed. That is we study the 
modification of the solution caused by incorporating 



e (A - *) md e - z ) m 

in (p0|) and fl3"2| ) respectively, just as the corresponding 9 functions were included in (|28| ) and 
In the small z, large Q 2 regime the kinematic constraint gf, < (1 — z)k\j z is a stronger 
limitation than the angular ordering constraint q 2 < Q 2 /z 2 , and we anticipate that the CCFM 
solution F(x, k?, Q 2 ) will become independent of Q 2 . In other words in this limit the kinematic 
constraint automatically embodies the angular ordering constraint || and since the former is 
independent of Q 2 the unintegrated gluon distribution F does not depend on this variable either. 
However, as Q 2 decreases below k\ the angular ordering constraint becomes stronger and F 
begins to decrease. We illustrate the effect in Fig. 4 which shows F versus Q 2 at different values 
of x for k\ = 10 GeV 2 . Here we have calculated F by numerically solving the CCFM equation 
( |30| ) as described in ref. [II], and then repeated the calculation with the kinematic constraints 



38]) incorporated. We also use these solutions to show in Fig. 5 the effective exponent A e ff of 
the integrated gluon distribution 



_ dln(xg) 

A e ff = 



<91n(l/x) 



rQ 2 

with xg{x,Q 2 ) = / Fdk 2 T (39) 



as a function of x for fixed values of Q 2 . The gluon distributions are generated from a flat input 
(_p(°) constant as x — > 0), and so the rapid rise with decreasing x(xg ~ x~ Acff ) is generated 
by perturbative QCD resummation effects via the CCFM equation. As for the BFKL equation 
with fixed as, we notice from Fig. 5 that for the CCFM equation A e fr becomes smaller (by 
about 0.1) on imposing the kinematic constraint. 



4. Impact on observables 

Predictions for observables at small x are driven by the behaviour of the gluon distribution, 
since the gluon is by far the dominant parton in this regime. Here we study three quantities, 
which are being measured at HERA, that are especially sensitive to the gluon distribution 
at small x; namely dF2/dlogQ 2 , the charm component of Fi and high energy diffractive J/ip 
photoproduction. The first two observables are linearly dependent on the gluon density, whereas 
the third has a quadratic dependence which considerably increases its sensitivity to the gluon. 

In deep inelastic scattering the virtual photon couples to the gluon via the g —>■ qq transition. 
We therefore calculate the structure function F2 from the unintegrated gluon distribution F 
using the /c^-factorization theorem 

F 2 (x,Q 2 ) = Y,Jdk 2 T ^ J d 2 KF(x',k 2 T ,K 2 ) F g box ^,K T ,k T ,Q 2 ,m q ^) + if (40) 
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where i^ box includes both the quark "box" and "crossed box" contributions which originate from 
virtual photon-virtual gluon qq production, that is from 7$ — > qq. The convolution is sketched 
in Fig. 6. For the u, d and s quark contributions we take the quark mass m q = 0, while for the 
charm component we take m c = 1.5 GeV. The explicit expressions for F box including quark 
mass effects can be foundf] in ref. fll3"l ; the argument of as in F^ ox is taken to be (k|> + m 2 ,) 
where Wq = 1 GeV 2 for u, d, s quarks and m 2 , = m 2 for the c quark. The results are not very 
sensitive to variations of mo about these values. 

The background contribution F 2 ~ Fzix, Q 2 ) at large x, but is a slowly varying function of 
x and Q 2 at small x. For example we extrapolate below x = 0.1 using the soft Pomeron x~ 008 
behaviour. However, predictions for the slope, dF 2 / d\nQ 2 are particularly insensitive to any 
ambiguities due to F 2 ■ The small x approximation of the CCFM equation, (|3"0"D, that we have 
used amounts to setting the Sudakov form factor As = 1 and to approximating the gluon-gluon 
splitting function by its singular term as z — > 0, that is P gg ~ 6/z. A 5 represents the virtual 
corrections which cancel the singularities at z — 1. To make a realistic comparison with data, 
we allow for the remaining finite terms in P gg by multiplying the solution F(x, k 2 ^, Q 2 ) by the 
factor 

/ 33 + 272/ fQ 2 _ . 2 dq 2 \ 

where the number of active flavours nj = 4. At large x(x > 0.1) and small Q 2 (Q 2 ^ 20 GeV 2 ) 
the CCFM solutions, with and without the kinematic constraint imposed, are found to agree 
with each other, and also with the solution obtained from the double-leading-logarithm (DLL) 
approximation of ( j30|) in which we replace 0(Q — zq) by Q(Q — q) and set A# = 1. However, 
at larger Q 2 (beyond the range of the data that we consider here) some care is needed. The 
kinematic constraint is only applicable in the small x region and so the normalisation of the 
gluon is suspect at large x, particularly for large Q 2 . We therefore renormalise the solution 
of the CCFM equation with the kinematic constraint imposed so as to agree for x > 0.1 
with the unmodified solution and its DLL approximation. In this way we allow for the small x 
approximation of the equation. The renormalisation only affects the solution for Q 2 > 20 GeV 2 . 

In Fig. 7 we show the predictions for F2(x,Q 2 ) at small x together with the latest HERA 
measurements. Including the kinematic constraint (p3|) in the CCFM equation for the gluon 
F(x, k 2 -,, Q 2 ) has the effect of taking us from the dashed to the continuous curves in Fig. 7. 
The relevant comparison is the slope (dF 2 / 'dlnQ 2 ) of the curves which is proportional to the 
gluon distribution. With the present experimental errors the comparison is inconclusive, but 
it is evident that, as the statistical and systematic errors are reduced, future measurements of 
<9i<2 /<9 In Q 2 will give insight into the properties of the gluon distribution F(x, k\, Q 2 ). 



Recently the charm component of F 2 has been measured |TJ] at HERA in the small x 
region. These measurements of F 2 (x, Q 2 ) are shown in Fig. 8, together with earlier EMC 
values [To"] at larger x. We compare these data with the charm component F^ determined from 

6 There is a typographical error in the expression for x' below eq.(19) in ref. [Q; the factor — f3) should 
be in the denominator. 
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the unintegrated gluon distribution, F, using the c quark contribution to the fc^-factorization 
formula (40). We show the values obtained by taking the mass of the charm quark to be 
m c = 1.4 and 1.7 GeV. We also show predictions based on GRV [I5| and MRS partons [ 17 . 
The first of these two is obtained from jg — > cc at NLO |Tj| using massive charm quarks and the 
integrated (GRV) gluon distribution. In the MRS analyses [|17], |l^] the charm quark is treated as 
a parton. The charm distribution is assumed to be zero for Q 2 < m 2 , while above this threshold 
(Q 2 > in 2 ) it is evolved assuming that m c = 0. The value m 2 = 2.7 GeV 2 is determined by 



fitting to the EMC data |L5| for F 2 C . Although the HI small x data are preliminary, it is clear 
that an improved measurement of i 7 ^ will be valuable. At present there are indications that 
our /^-factorization approach underestimates the HI data at the lower Q 2 values; in fact the 
imposition of the kinematic constraint worsens our previous description of these data |20| . 



The third observable process that we study is high energy diffractive J /if) photoproduction, 
7P — > J/ ipp. It offers a sensitive probe of the gluon at small x JTL], |22[. The amplitude can 
be factored into the product of the 7 — > cc transition, the scattering of the cc quark pair on 
the proton via (colourless) two-gluon exchange, and finally the formation of the J /if) from the 
outgoing cc pair. The crucial observation is that at high energy the scattering on the proton 
occurs over a much shorter timescale than the 7 — > cc fluctuation or the J/ip formation times. 

Since diffractive J /if) photoproduction is essentially an elastic process, the cross section is 
dependent on the square of the unintegrated gluon distribution, F(x,k^,Q 2 ). The relevant 
values of x and Q are x = Ml/W 2 and Q = AfJ/4, where Afy is the mass of the J /if) meson 
and W is the ^p centre-of-mass energy. The cross section is given by |23 

1 da 



ofyp-* J/ipp) 



b dt 



vr 3 MX(Q 2 )r e 



36a 



M 2 + Ak^ / 



F(x, k 2 , 



(41) 



where T ee is the leptonic width describing the J /if) — > e + e~ decay, a is the QED coupling, 
and b is the slope parameter of the differential cross section, da/dt = Aexp(— b\t\). We take 
the experimental value b = 4.5 GeV -2 . We include the effects of cc rescattering and the small 
contribution of the real part of the amplitude as described in ref. [^3j. It was noted in ref. [p3|| 
that the effects of Fermi motion of the c and c quarks in the J /if) lead to a sizeable (±30%) 
uncertainty in the normalization of the perturbative QCD prediction of the photoproduction 
cross section, but that the "shape" of the W (or x) dependence is unaffected. 

The predictions for diffractive J /if) photoproduction are compared with recent HERA data 
in Fig. 9. At present the data extend up to energy W ~ 140 GeV, that is down to x ~ 5 x 10 -4 . 
The prediction in the absence of the kinematic constraint (the dashed curve) implies that the 
gluon increases too fast with decreasing x. On the other hand if the kinematic constraint is 
incorporated in the CCFM equation, then the continuous curve is obtained and the description 
is improved. 
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For completeness we also show in Fig. 9 the description of the J/ ip data calculated from two 
recent sets of partons (GRV [T6| and MRS(A') |T7| ) as described in ref. j2"3fl . Neither parton set 
incorporates ln(l/x) resummation effects. The J/ip data appear to favour the phenomenological 
gluon distribution of the latter set of partons. 

5. Conclusions 

We have studied the behaviour of the gluon distribution of the proton in the small x regime, 
10~ 4 < x < 10~ 3 , that has recently become accessible to the experiments being performed at 
the HERA electron-proton collider. In this regime it is necessary to work in terms of the unin- 
tegrated gluon distribution F(x, k%) and to resum ln(l/a;) contributions. To leading order, the 
resummation is accomplished by the BFKL equation. The a^ ln n (l/:r) contribution corresponds 
to an effective n-rung ladder diagram arising from a space-like chain of n gluon emissions. The 
solution F(x, h?) of the BFKL equation shows, with decreasing x, two characteristic features. 
First an x~ x growth (where A = 0:54 In 2 for fixed as, or A ~ 0.5 if as is allowed to run), and 
second, a diffusion or random walk in In k\ as we proceed along the gluon chain. 

The full next-to-leading order summation of terms ln ra_1 (l/:r) is not yet known, but we 
would expect the description of F 2 to be sensitive to this correction. However, two important 
higher-order effects can already be investigated. One effect is due to the angular ordering of 
gluon emissions, which leads to the CCFM, rather than the BFKL, equation for the gluon 
distribution. The solution now depends on an additional scale that is required to specify the 
maximum angle of gluon emission, which turns out to be essentially the scale Q 2 of the probe. 
That the solution F(x, k 2 -*, Q 2 ) will be dependent on Q 2 is indeed evident from the angular 
ordering constraint Q(Q — qz) in the CCFM equation, (|30|). At very small x the constraint is 
automatically satisfied and the CCFM solution reduces to the BFKL form F(x, kj). Numerical 
solutions to the CCFM equation were obtained in ref. [T1J . 

The second higher-order effect is due to the imposition of the kinematic constraint Q(k^/q^— 
z) which is required for the validity of the BFKL or CCFM equation at small x. The constraint is 
needed to ensure that the virtuality of the gluons along the chain is controlled by the transverse 
momenta, that is \k 2 \ ~ k\. The major aim of this paper is to explore the consequences of 
implementing this constraint. For fixed as the introduction of the constraint preserves the 
scale invariance of the BFKL equation. The x~ x behaviour of the solution as x —>■ can 
therefore be obtained using analytic methods. We found that the BFKL intercept A = as A In 2 
is significantly reduced. The details are shown in Fig. 2. We see that it is insufficient to consider 
just next-to-leading effects; higher-order effects are important. For running we solved the 
CCFM equation numerically and obtained the gluon distribution F(x, k\, Q 2 ) with and without 
the kinematic constraint imposed. In this case the constraint has the effect of reducing A by 
about 0.1, see Fig. 5. For Q 2 > k\ the kinematic constraint is more severe than that due to 
angular ordering and so the gluon distribution becomes independent of Q 2 , see Fig. 4. 

In Section 4 we studied the impact of imposing the kinematic constraint on the description 
of three observables which are sensitive to the gluon distribution at small x and which are 
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being measured at HERA. The observables are dF 2 /d\nQ 2 , F£(x, Q 2 ) and the W dependence 
of the cross section for high energy diffractive J/ip photoproduction. The effects are significant 
in particular regions of phase space; contrast the dashed and continuous curves in Figs. 7 
and 9. As expected J/ip photoproduction offers an especially sensitive measure of the gluon. 
These comparisons with data should be regarded as exploratory, in other words, our calculation 
should be viewed as a preliminary step towards a full theoretical treatment of higher-order 
ln(l/x) contributions. However, it is clear that the effects discussed in this paper should be 
incorporated in any realistic analysis of improved small x data. 
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Figure Captions 



Fig. 1 Gluon emission which forms the basis of the BFKL equation (|T]) for the unintegrated 
gluon distribution F(x,kj>)- x and x/z are the longitudinal momentum fractions of the 
proton's momentum carried by the respective gluons. Throughout we use qx and k^ to 
denote, respectively, the transverse momentum of an emitted gluon and of a gluon along 
the chain. 

Fig. 2 The exponent A of the x~ x behaviour of the gluon distribution obtained by solving the 
BFKL equation (a) with (continuous curve) and (b) without (dashed curve) the kinematic 
constraint imposed, as a function of (fixed) as = 3as/n- The dashed curve is A = as4 In 2. 
The dotted curve (c) is the value of the exponent that is obtained if we keep only the 
next-to- leading order modification due to the kinematic constraint. 

Fig. 3 The continuous curves are the effective exponent A e fj = d\n.F/d\xi{l/x) calculated from 
the numerical solution F(x, kj) of the BFKL equation which incorporates the kinematic 
constraint, for k\ = 1000, 100, 10 and 4 GeV 2 . The coupling a s = 3a s /n = 0.2. The 
dashed lines indicate the value of the exponent of the x — > analytic solutions of the 
BFKL equation with and without the kinematic constraint included, that is the values of 
curves (b) and (a) of Fig. 2 respectively at a~s = 0.2. 

Fig. 4 The Q 2 dependence of the gluon distribution F(x, kj-, Q 2 ) obtained by solving the CCFM 
equation with (continuous curves) and without (dashed curves) the kinematic constraint 
(1331) included. Results are shown for x = 10~ 5 , 10~ 4 , 10~ 3 and 10~ 2 . 



Fig. 5 The effective exponent A c ff = d\a{xg)/d\sx{l/x) obtained by solving the CCFM equa- 
tion with (continuous curves) and without (dashed curves) the kinematic constraint (|33| ) 
included. Results are shown for Q 2 = 1000, 100, 10 and 4 GeV 2 . 

Fig. 6 Pictorial representation of the kr factorization formula, that is of the convolution F% = 
J2 q F (g> F 9 box of (pE0"P. F{x', fey, k 2 ) is the unintegrated gluon distribution and Y, q Fq OX is 



the off-shell gluon structure function, which at lowest order is determined by the quark 
box (and "crossed box") contributions. 

Fig. 7 Predictions of the proton structure function 1*2(2, Q 2 ) as a function of InQ 2 , at fixed 
values of x, compared with recent measurements made by the experiments at HERA 
|1], 0] . The continuous and dashed curves are respectively the predictions obtained from 
the CCFM equation, via the /^-factorization theorem, with and without the kinematic 
constraint ( |33| ) incorporated. The dotted curves are the predictions obtained from the 
GRV set of partons |T6[ . 



Fig. 8 Predictions for the charm component F 2 C of the proton structure function, F 2 , compared to 
recent preliminary HI measurements [l4j] and older EMC data |TJ[ . The predictions were 
obtained by solving the CCFM equation (with kinematic constraint) for the unintegrated 
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gluon and then using the /^-factorization formula with m c = 1.4 GeV (upper continuous 
curve) and m c = 1.7 GeV (lower continuous curve) Also shown are the next-to-leading 
order predictions [T^| based on GRV [|TjJ] and MRS [T7| partons. 



9 The measurements ]3| of the cross section for diffractive J/ip photoproduction compared 
with the perturbative QCD description based on the unintegrated gluon distribution 
obtained by solving the CCFM equation with (continuous curve) and without (dashed 
curve) the kinematic constraint fl33|) included. The dotted and dash-dotted curves are 



the predictions obtained from the GRV and MRS (A') set of partons [16|, [L7]], calculated 
as in ref. [E3|. 
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